Topology and manipulation of multiferroic hybrid domains in M11WO4 
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An investigation of the spatially resolved distribution of domains in the multiferroic phase of 
MnW04 reveals that characteristic features of magnetic and ferroelectric domains are inseparably 
entangled. Consequently, the concept of multiferroic hybrid domains is introduced for compounds 
in which ferroelectricity is induced by magnetic order. The three-dimensional structure of the 
domains is resolved. Annealing cycles reveal a topological memory effect that goes beyond previously 
reported memory effects and allows one to reconstruct the entire multiferroic multidomain structure 
subsequent to quenching it. 



In the field of strongly correlated electron systems ma- 
terials with cross-correlated magnetic and electric prop- 
erties, called magnetoelectrics, are intensely discussed be- 
cause of their potential for controlling magnetic prop- 
erties by an electric voltage. Among them, the mag- 
netoelectric multiferroics may be most prominent be- 
cause due to a coexistence of magnetic and electric or- 
der they can developparticularly pronounced magneto- 
electric interactions. [H, [i[ Recently it was demonstrated 
that intrinsically strong ("giant") magnetoelectric ef- 
fects are present in the so-called joint-order-parameter 
multiferroics in which magnetic long-range order breaks 
the inversion symmetry and induces a spontaneous 
polarization. 0,11,01 For example, in TbMn0 3 , Ni 3 V 2 8 , 
and MnW04, a spiral arrangement of spins violates the 
inversion symmetry and causes a spontaneous polariza- 
tion 

P oc eij x (Si x Sj) , (1) 

with as unit vector connecting neighboring spins at 
sites i and j and (Si x Sj) as vector chirality. 0, 0, 
The magnitude and direction of P are determined by the 
magnetic order only, so that a unique correlation between 
the magnetic and ferroelectric (FEL) order parameters is 
obtained. Although the spontaneous polarization is usu- 
ally small, its robustness [9] renders joint-order-parameter 
multiferroics interesting for future applications. 

An essential feature of any ferroic material is the pres- 
ence of domains. They determine the switching of in- 
formation bits in memory devices and the technologi- 
cal performance of permanent magnets. At its root, 
any magnetoelectric interaction in a multiferroic corre- 
sponds to an interaction of its magnetic and electric do- 
mains. Hence, understanding giant magnetoelectric ef- 
fects means understanding the nature and interactions 
of the multidomain state in the joint-order-parameter 
multiferroics. [16] However, although domains are known 
to be present in these compounds by a variety of (mostly 
indirect) experiments 0, El, EH El, El, El, EH the prime 
target of previous investigations was actually to remove 



these domains by converting the samples to or in between 
single-domain states. 

Therefore, none of the existing publications addresses 
the three-dimensional distribution of domains in a joint- 
order-parameter multiferroic and the effects guiding their 
formation. With this, essential questions regarding the 
nature of the domains remain unclear. For example, to 
what extent can a domain actually be called ferroelectric 
if the spontaneous polarization is magnetically induced? 

In this report, the three-dimensional topology of do- 
mains in the joint-order-parameter multiferroic MnWCU 
is resolved by optical second harmonic generation (SHG). 
The spatial distribution of the domains, their response to 
external fields, and annealing procedures reveal some fea- 
tures that are uniquely associated to a magnetic domain 
state and others that point unambiguously to ferroelec- 
tric domains. The concept of "multiferroic hybrid do- 
mains" is thus introduced whereas a description in terms 
of "magnetic domains" or "electric domains" or even 
"magnetic domains coexisting with electric domains" is 
no longer appropriate. 

The evolution of the multiferroic phases and their 
domains in joint-order-parameter compounds can be 
described by Landau theory. In the case of spin- 
spiral systems it was shown that a coexistence of two 
magnetic order parameters is required for obtaining 
multiferroicity.0, El, El, The first magnetic order 
parameter generally guides the system into a phase with 
an incommensurate non-polar spin arrangement. Upon 
cooling another magnetic transition described by a sec- 
ond magnetic order parameter evolves, so that the com- 
bination of both breaks the inversion symmetry of the 
crystal. As a consequence, a spontaneous polarization 
according to Eq. (pQ) can emerge. 

For the investigation of the local structure of domains 
in spin-spiral ferroelectrics the choice of MnWCU as a 
model system suggests itself. On the one hand, its mag- 
netic lattice is rather simple, because no rare-earth and 
just one kind of transition- metal ion contributes to the 
magnetic order. On the other hand, the rich magnetic 
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FIG. 1: (a) Magnetic phase transitions in MnWCU. Below the 
Neel temperature Tn the Mn 3+ moments display collinear in- 
commensurate long-range order within the easy plane. At T2 
an additional transverse spin component orders which leads 
to a helical incommensurate arrangement of spins. In the 
magnetic ground state below T\ a commensurate up-up-down- 
down spin ordering is obtained. The spin spiral in the AF2 
phase induces a spontaneous electric polarization along the 
y axis while the AF1 and AF3 phases are not multiferroic. 
(b - d) Spatially resolved SHG measurements in the AF1 to 
AF3 phases. A domain structure (with black lines reveal- 
ing domain walls) is detected in the multiferroic AF2 phase 
only, x, y, and z represent a Cartesian coordinate system, ap- 
proximating the monoclinic unit cell with lattice parameters 
a = 4.83 A, b = 5.76 A, c = 4.99 A and f3 = 91.1° « 90°. 



phase diagram includes transitions from the multiferroic 
phase to neighboring ordered or disordered states which 
allows one to probe memory effects in either. 0, [2l| 

Three magnetically ordered phases determine the low 
temperature behavior of MnWC^. In the AF3 phase 
right below Tn = 13.5 K incommensurate antiferromag- 
netic (AFM) ordering of the Mn 3+ moments described by 
the order parameter (9af3 occurs. The spins align in a 
collinear way, pointing along the easy-axis within the xz 
plane so that they enclose an angle of about 34° with the 
x axis of the monoclinic crystal as illustrated in Fig.QJa). 
Their magnitude is sinusoidally modulated, leading to a 
two-dimensional incommensurate spin-density wave with 
k= (-0.214, \, 0.457). 

In the multiferroic AF2 phase below T2 = 12.7 K an 
additional transverse spin component orders and the spin 
density wave becomes an elliptical spin spiral while re- 
taining k = (—0.214, |, 0.457) as wave vector. The re- 
lated magnetic order parameter 0af2 evolves through 
a second-order phase transition at T2 and coexists with 
0AF3- As detailed above, simultaneous presence of the 
two magnetic order parameters breaks the inversion sym- 
metry and induces a spontaneous polarization, here 



P y OC (9AF30AF2 



(2) 



along to the y axis, and establishes multiferroicity in 
MnW04. Note that Oaf2 can be reoriented below T2 
while Oaf3 is frozen. [H, [H, Thus, any reversal of 
&AF2 is coupled one-to-one to a reversal of P v . [23L l24j 

In the AF1 phase, MnWCU displays collinear com- 
mensurate AFM order. The associated wave vector is 



k = (±~, ~, |) and describes an up- up- down- down spin 
structure with spins aligned along the easy- axis. In 
agreement with Eq. (pQ), the magnetic order parame- 
ter 0AF2 and, thus, P y , are quenched at the first-order 
AF2 — > AF1 transition at T x = 7.6 K while G A F3 re- 
mains. 

Because of the intricate interplay of the order parame- 
ters 0AF2 and 0AF3 the number and types of domains for 
the three phases will differ. A powerful technique for in- 
vestigating domain structures in systems with magnetic 
or electric order is second harmonic generation (SHG). 
Optical SHG describes the induction of a light wave at 
frequency 2uj by a light wave at frequency uo. The pro- 
cess picks up the symmetry changes imposed by the long- 
range order by coupling directly to the corresponding or- 
der parameter and its different orientation in different 
domains. A detailed discussion of the technical aspects 
of SHG in ferroic systems in general [25| and in MnWCU 
in particular [l2| was already published so that we restrict 
ourselves here to the application of SHG for imaging the 
domain structure of MnWC^. 

Obviously, SHG contributions arising at Tn probe the 
magnetic order parameter Oaf3, while those emerging at 
T2 involve coupling to Oaf2 and, because of Eq. ([2|), to 
P y . As overview, Figs. [2(b)- (d) show spatially resolved 
measurements of the SHG intensity in the AF1, AF2, 
and AF3 phase, respectively. While the AF1 and the 
AF3 phases reveal a homogeneous distribution of SHG 
intensity, the multiferroic AF2 phase exhibits about ten 
curved black lines on an otherwise homogeneous back- 
ground. Such lines are a hallmark for the presence of 
domains with opposite orientation of the corresponding 
order parameter. Because of the linear coupling to the or- 
der parameter, the SHG light field experiences a sign re- 
versal corresponding to a 180° phase shift when crossing 
the domain wall. This leads to destructive interference 
in the vicinity of the domain wall and, therefore, to the 
black lines. According to Fig. [2 such 180° domains are 
present in the multiferroic AF2 phase only while being 
absent in the AF1 and AF3 phases. 

In order to understand this observation we have to an- 
alyze the symmetry of the AF1 to AF3 phases. The mon- 
oclinic paraelectric and paramagnetic phase of MnW04 
belongs to the point-group 2 y /m y V. This point symme- 
try is preserved by the commensurate up- up- down- down 
spin configuration in the AF1 phase and also by the in- 
commensurate spin-density wave in the AF3 phase be- 
cause the magnetic order breaks translation symmetries 
only. [22] Therefore, a formation of domains with differ- 
ent orientation of the order parameter should indeed not 
occur. HH In contrast, the multiferroic AF2 phase pos- 
sesses the point symmetry 2 y Y with half the number of 
point-symmetry operations as in the group 2 y /m y V . As 
observed, this leads to two domains with an opposite sign 
of the magnetic order parameter Oaf2 corresponding to 
an opposite magnetic vector chirality and, via Eq. ([2]), to 
an opposite sign of the spontaneous polarization P y . 

In the following, we will focus on the discussion of the 
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FIG. 2: Spatially resolved SHG images reveal the hybrid na- 
ture of the multiferroic domains in MnWCU- (a) Domain 
structure in the xz plane. A pronounced elongation along the 
magnetic easy axis of the crystal (yellow arrow) and the bub- 
ble topology reflect the magnetic aspect of the domains, (b, 
c) Domain structure in electric fields of (b) 167 kV/cm and 
(c) 217 kV/cm applied along the y axis, (d, e) Domain struc- 
ture in magnetic fields of (d) T and (e) 5 T applied along 
the x axis. Only the electric field affects the distribution of 
the domains, thus reflecting their electric aspect. 



180° domains in the multiferroic AF2 phase. First of all, 
we will investigate the topology of the AF2 domains and 
their response to external poling fields. Figure[2](a) shows 
the distribution of the AF2 domains in the xz plane after 
zero-field cooling. We find a labyrinth-like arrangement 
of dark and bright areas, being distributed with approx- 
imately equal proportions across the sample. The differ- 
ent level of brightness corresponds to opposite magnetic 
vector chirality and spontaneous polarization according 
to Eq. (pQ). In contrast to Fig.[T](c) a different brightness 



between opposite domains (in addition to the mere do- 
main walls) is observed because of the interference of the 
SHG light field arising from the multiferroic order with 
a homogeneous SHG light field generated by the crys- 
tal lattice. [25(| The topology of the domains in Fig. [TJa) 
reveals two preferential directions: along the z axis and 
along a line including an angle of about 34° with the x 
axis. 

Figure [21(a) clearly expresses the magnetic origin of 
the domain structure. First, the texture of the domains 
strikingly resembles the patterns of bubble and stripe do- 
mains, universally attributed to modulated phases with 
a preselected equilibrium periodicity. (2?J In MnW04 it 
is the long-range magnetic interaction that stabilizes the 
modulated AF2 phase and the periodicity is determined 
by the magnetic wave vector k. Second, the preferen- 
tial direction of the domain walls parallel to the arrow 
in Fig. [21(a) is in striking coincidence with the direction 
of the magnetic easy-axis of MnWGU which encloses an 
angle of 34° — 37° with the x axis in the xz plane. [22| 

In Figs.[2](b - e) the response of the domain structure in 
the yz plane to electric and magnetic fields is shown. It is 
obvious that with an electric field applied along the y axis 
small changes in the order of 1 kV/cm lead to pronounced 
changes in the topology of the AF2 domains. In contrast, 
a magnetic field of up to 5 T applied along the x axis does 
not alter the domain structure at all [Fig. [2Jd) and [21(e)], 
and the same holds for fields applied along the y or z 
axis. Figures [2Kb - e) thus emphasize the electric nature 
of the AF2 domain structure. 

We therefore conclude that the AF2 domains exhibits 
hallmarks of both a magnetically and an electrically or- 
dered state. The topology is determined by the magnetic 
character of the domains, whereas the field response re- 
flects the electric character. Because of this bilateral na- 
ture, only a denomination of the AF2 domains as multi- 
ferroic hybrid domains seems appropriate whereas a de- 
scription in terms of ferroelectric or antiferromagnetic 
remains incomplete. This is in stark contrast to ear- 
lier work[l2| where magnetic and ferroelectric domains 
where still considered as strictly separate entities. It was 
already pointed out that the improper nature of the fer- 
roelectric polarization implies rigid coupling to the mag- 
netic order parameters. 0, [Hf However, this seems to 
be the first work where the consequences of this relation 
on the level of the domain structure is revealed. Only 
this leads to the realization that an extended concept of 
multiferroic hybrid domains is required. 

The full three-dimensional topology of the multiferroic 
domains is presented in Fig. [3] for three differently ori- 
ented MnW04 samples obtained from the same batch. 
Figure [3](a) shows the domain structure in the xz plane. 
Like in Fig. [21(a), dark and bright regions correspond to 
the two possible multiferroic domains with opposite order 
parameters. We can see the aforementioned elongation 
of domains along the magnetic easy-axis with the sponta- 
neous polarization P y pointing into and out of the plane. 
The propagation of domain walls along the direction of 




FIG. 3: Three-dimensional distribution of the multiferroic do- 
mains, (a - c) Domains in the xz, yz, and xy plane of MnWCU 
samples taken from the same batch, (d) Three-dimensional 
visualization of the multiferroic domain structure in (a - c). 



the spontaneous polarization is revealed by Figs. [31(b) 
and [31(c). Interestingly, the domain walls, here again in- 
dicated by the black lines, continue rather straight along 
the y axis of the crystal. The domains tend to form 
platelets in planes defined by the magnetic easy-axis and 
the direction of the spontaneous polarization. With a lat- 
eral extension evaluated as 170 ± 60 /im[28| the domains 
are surprisingly large. For an improved visualization of 
the anisotropic multiferroic domain structure, Fig. [3](d) 
shows a three-dimensional simulation based on the dis- 
tribution of domains in Figs. [3](a - c). The respective 
domain structure were projected onto three faces of a 
cuboid and mended minimally at the edges. 

In the following we consider the dynamic aspects in 
the formation of the domains by applying annealing pro- 
cedures across the boundaries limiting the multiferroic 
phase. First of all, we consider the effect of temperature- 
annealing in an order <-> disorder cycle from the AF2 
phase to the paramagnetic state above T/v and back to 
the multiferroic phase. Initially, the sample was zero- 
field-cooled from room temperature into the AF2 phase 
which leads to the SHG image shown in Fig. Ufa). In 
agreement with Fig. HJ a multiplicity of domains is ob- 
tained. Subsequently, the order <-> disorder annealing 




FIG. 4: Response of the multiferroic domains to a tempera- 
ture annealing cycle through the paramagnetic phase, (a) Do- 
main structure after initial zero-field cooling from room tem- 
perature to the multiferroic AF2 phase, (b) Domain structure 
after subsequent application of the annealing cycle. 



cycle was applied which leads to the SHG image shown 
in Fig. Hfb). Most of the domain walls have vanished and 
only a small fraction of the MnW04 remains in a domain 
state opposite to the rest of the sample. Apparently, the 
annealing procedure tends to drive the sample towards 
a single domain state which was confirmed by repeating 
the experiment more than ten times. 

Figure [H leads to two conclusions. First, no mem- 
ory effect is observed in the order <-> disorder anneal- 
ing cycle. Second, while conventional ferroelectrics tend 
to form a multiplicity of domains for minimizing electric 
stray fields, an ideal antiferromagnets tends to approach 
a single-domain state. Hence, just like the domain 
structure itself (Fig. [2]) the dynamic domain topology is 
dominated by the magnetic aspect of the hybrid multi- 
ferroic order. 

In the second step, we investigated the effect of an or- 
der <-> order annealing cycle from the AF2 phase to the 
AF1 phase and back to the multiferroic state. This tran- 
sition is of particular interest because of a FEL memory 
effect reported earlier: [2l|, [30[ A single-domain state in 
the multiferroic phase is memorized in the non-polar AF1 
phase and reemerges when reentering the AF2 phase. 
However, all current data were gained by integral tech- 
niques such as pyroelectric current measurements, so that 
it is not known to what extend the memory effect applies 
to the domain structure of the multiferroic multidomain 
state. 

This is investigated in Fig. [5] which compares the do- 
main structure of the AF2 phase before [panels (b), (d)] 
and after [panels (c), (e)] the annealing cycle through the 
AF1 phase. Figure 03(a) illustrates that the phase bound- 
ary between the AF2 and the AF1 phase can be crossed 
by magnetic-field or temperature tuning. Figures [51(b) 
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FIG. 5: Response of the multiferroic domains to magnetic- 
field (H) and temperature (T) annealing cycles through the 
AF1 phase, (a) Sketch of the (if, T) phase diagram with ar- 
rows indicating the respective annealing procedure, (b, d) 
Domain structure after initial zero-field cooling from room 
temperature to the multiferroic AF2 phase, (c, e) Domain 
structure after application of the respective annealing cycle, 
(f) Temperature dependence of the SHG signal (xyxx compo- 
nent at 1.95 eV)[T2] in a temperature decreasing run. Nonzero 
SHG yield in the AF1 phase points to a residual polar contri- 
bution in the AF1 phase as probable basis of the topological 
memory effect revealed in panels (b - e). 



and [He) reveal the effects of thermal annealing. Minor 
changes in the roughness of the domain walls and, for a 
small fraction of the walls, shifts in the order of 100 /am 
are observed but the general distribution of the domains 
is preserved. Annealing in a magnetic field along the y 
axis reveals no changes at all in the domain structure. 
Hence, not only a single-domain state is memorized in 
the non-chiral non-polar AF1 phase — actually the en- 
tire topology of a multidomain state is preserved. The 
memory effect is thus much more rigid than established 
up to now. 

Figure [5] immediately raises the question for the origin 
of such a memory effect. Pinning of the domain structure 
by structural defects can be excluded. This mechanism 
would also preserve the domain structure in the order <-» 
disorder transition to the paramagnetic phase contrary 
to what Fig. [4] shows. For revealing its origin we mea- 
sured the temperature dependence of the magnetic order 
parameter Oaf2 across the AF2 — > AF1 transition. As 



shown in Fig 03(f), the SHG contribution that is related 
to the magnetic order parameter Oaf2 does not vanish 
abruptly at the phase boundary to the AF1 state. In- 
stead, it diminishes gradually with temperature after an 
initial step-like decrease. The remanent SHG signal re- 
veals a coexistence of the AF1 and AF2 phase within a 
broad temperature region, enabled by the first-order na- 
ture of the transition. Hence, residual nuclei of the AF2 
phase explain the "polar memory" of MnWO4.[30f The 
distribution of the nuclei must be homogeneous and dense 
in order to explain pinning of the entire topology of a mul- 
tidomain state. Note that 03(f) is the first measurement 
directly proving the existence of polar inclusions within 
the AF1 phase, which is possible because of the extraor- 
dinary sensitivity of SHG to nanoscopic inclusions. [3lj 

In summary, we revealed that the magnetically induced 
ferroelectric phase of the joint-order-parameter multifer- 
roic MnW04 forms novel types of domains for which 
characteristic features of magnetic and ferroelectric do- 
mains are inseparably entangled. Here the denomination 
as multiferroic hybrid domain is introduced. In MnW04 
the topology of these chimera domains is determined by 
their magnetic nature, whereas the field response reflects 
the electric character. The three-dimensional distribu- 
tion of the domains was investigated and annealing cy- 
cles revealed a topological memory effect allowing one to 
reconstruct the entire multidomain structure subsequent 
to quenching it. 

The present work should provide a basis for the devel- 
opment of a model explaining the topology of domains 
in compounds with magnetically induced ferroelectric- 
ity. Regarding long-term application, a vastly increas- 
ing variety of joint-order-parameter multiferroics allow- 
ing ferroelectricity according to Eq. ([2j) are at our dis- 
posal for multiferroic domain control. The discovery of 
mechanisms promoting magnetically induced ferroelec- 
tricity different from Eq. ([2]) is only a question of time 
and may eventually lead us towards high-temperature 
applications. Because of the bilateral nature of the mul- 
tiferroic hybrid domains systems with conical spin spi- 
rals allow one to exert rigid electric-field control of a 
macroscopic magnetization. [23j Since magnetization and 
polarization are manifestations of the same multiferroic 
hybrid domain state the electric field will always act 
on both. Rapid magnetization reversal with ultrashort 
electric-field pulses may thus become feasible. 
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